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ABSTRACT 

A computer based study to assess the feasibility of agricultural 
utilization of reject heat from a nominal 1000 MWe power plant began in 
September 1973 ly an interdepartmental team at MSU. This paper reviews 
and interprets the results of this systems study. Components used 
include numerical models representing three agricultural subsystems : 
pond heating for catfish culture, field warming for vegetable production, 
and greenhouse heating for ornamental flowers as well as models for 
the power plant thermal discharge and a cooling reservoir for the 
dissipation of heat not usable by the agricultural subsystems. The 
components of the system were integrated by optimizing a net present 
value criterion. Ownership-management options for operating the 
agriciiltural subsystems are siiggested. This paper presents impressions 
on the usefulness of the factors included in the methodology \ised. 

Future research needs for this methodology are outlined. 


INTRODUCTION 

Over the last decade there have been numerous conferences to assess the 
potential for and problems with utilizing reject heat in beneficial 
ways. The early conferences usually had two separate components; (1) 
a collection of technical papers with each paper on a specific aspect 
of waste heat utilization, and (2) a concluding statement or set of 
recommendations reiterating the diverse factors that need to be 
considered in a unified assessment. When discussions were held at 
Michigan State University (MSU) in September 1973 to institute a 
feasibility study on waste heat utilization in Michigan, the consensus 
of opinion was that (1) adequate technical data for Michigan conditions 
could be developed from the results of field experiments conducted 
elsewhere, and (2) that a computer based approach could be used to get 
a low cost, rapid, overall assessment. The key feature of the proposed 
study was the interaction over the duration of the study by the small 
but multidisciplinauy group. The group would be forced — by the necessity 
of developing mathematical models and analysis techniques — to isolate 
the dominant factors in previous research. 

Preceding page blank 


IPS 



VI-A-4 


The tentative orgajilzation of the study was defined in the following 
way. The candidate agricultural uses of waste heat were selected 
from systems studied "by others [.^,2] with giiidance from Dr. Laxry 
Boersma Cl] including preprint material from his research then in 
progress . Adequate data seemed to be available to model pond 
warming for freshwater fish culture, greenhouse heating, soilwarming 
for vegetable production, and several uses. Because the initial 
goal in the computer study was to develop an analysis "tool", the 
study was limited to these three agricultural uses. The work by 
Boersma *s group l_4,p.243] showed that important parameters in 
several components of a combined waste heat system coiald be defined, 
but they did not have time available to develop an integrated 
system. A combined system is one consisting of several components 
that are Independent. An integrated system is one where the size 
of the components, or their arrangement, or their operation is 
chosen to optimize a criterion. Previously Price had integrated 

the sizes of three components (fish pond, greenhouse, and recreational 
laJce) in a system using the criterion of maximum overall size of 
the component uses. It was optimized by a sensitivity study involving 
heat dissipation models and actual weather data. It was suspected 
that the cost of transporting the enoimious quantity of low temperature 
water to the components of the agricultural system might be ein 
important parameter Also a detailed economic analysis of 

the cost and revenue in the agricultural components would certainly 
be Important ‘to agribusiness \_9], Finally, a net present value 
criterion was taiken as the criterion that would be used to 
integrate the three agricultural subsystems. 

The eleven university personnel involved with the project for varying 
periods of time presented their preliminary conclusions based on 
this criterion in December 1975 During the first part of 1976, 

two of the personnel continued to made refinements on the study [ll]* 

Examples could be given of papers in ther literature where waste heat 
computer models were used for the following goals : (1 ) development of 
mathematical models that accurately represent the collected experimental 
data for specific agricultural uses, (2) development of detailed 
design specifications for a specific subsystem, (3) approximate 
evaluation of an informative criterion, and (4) computational 
proced\rres to optimize some criterion. The essence of the MSU stiidy 
involves the last two items. Informative information includes; 
cooling system net present value, enterprise shadow prices, agri- 
business net present value, and environmental concerns. 

Separable programming and a net present value criterion were a 
useful way to study the agriciILtural utilization of reject heat. 
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NET PRESENT VALUE 

Net present value (NPV) is used by the electric utility industry to 
compare alternative conventional cooling systems. It is used to justify 
decisions on the cooling system to their investors aind to state 
regulatory agencies. 


The net present value method accounts for the time value of money 
in combining the investment in initial capital and annual costs. 
The NPV was used as the criterion that was optimized to define the 
size and temperatures of the agricultiaral subsystems. For the 
agricultural subsystems the net present value is of the form; 
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Net present value of the integrated subsystems, positive 
terms j for properly chosen subsystems, 

Number- of subsystems, agricultimral plus reservoir. 

Area or size of subsystem j , acres 
Annual gross revenue for subsystem j. 

Annual cost for subsystem j. 

Discount rate or opportunity cost for capital. 

Initial capital outlay for subsystem j. 

Life of project or planning horizon, years, 

Index for time, 

Index for subsystem. 


While this objective function is usable for determining the optimal 
subsystem sizes, it does not allow comparisons with conventional 
cooling methods because it does not accoiant for the cost in transporting 
the warm water to the subsystems. The NPV for the water transport 
system is of the form; 
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additional variables are 

= Net present value of water transport system, or 
of conventinnal cooling system, always negative. 

= Capital outlay for water transport system, 

= Maikeup power and other annual costs for water transport. 


Ihe same form for Z* is used to compute the net present value of 
conventional cooling, with different suitable values for P* and K* 
of covirse. For the integrated aigricultural system, Z = Z* + Z**. 
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Ihe optimal integrated system contained 100 acres of tomato fields, 

160 acres of catfish ponds, a 350 acre cooling ledce, and no green- 
houses. Bakker et al. C^Oj concluded that the NPV for this integrated 
aigricultural system makes it competitive with conventional cooling 
methods. Ihe slightly revised NPV comparisons from Meekhof are 

shown in table 1 . The NPV for the integrated agricultviral system is 
less negative than the corresponding value for conventional cooling. 
Although the NPV is positive for each agriciiltural system, the value 
for the cooling reservoir has a slightly greater magnitude and is 
negative. When the NPV for the water tiansport shown in table 6 is 
included the total NPV is about as negative as that for a conventional 
cooling tower. 


SHADOW PRICES 

One of the programming techniques used to optimize the integrated 
system was linear programming (LP). Ihe IP output gives shadow 
price information useful to agribusiness in decision making. 

Shadow prices for activities (subsystems) indicates the effect on 
NPV of expanding the size of an enterprise. Shadow prices for 
constraints indicate the effect on the NPV for all subsystems of 
removing marketing and other constraints. Positive values are an 
indication of what entrepeneurs would pay for one additional unit. 
Table 3 [from ll] gives Scimple shadow prices for the optimal 
system. The negative value for waste heat in September indicates 
that adding a subsystem operating in that month would be desirable. 
Ihe magnitudes of the other sheidow prices show that the fish 
enterprise reso\irce is the major limiting factor. An explicit 
constraint in the LP program limited the catfish enterprise to 
less than 160 acres. 

Several resovirces are limited in the waste heat system and are 
represented by explicit and implicit constraints in the LP problem. 
Available markets limit the sizes of the tomato and catfish 
enterprises. The total quantity of waste heat is limited by the 
size of the electric power plant, assumed to be 1000 MWe. Because 
the reject heat is used to keep the subsystem temperatures above 
specified levels, the heat dissipation considerations discussed 
later imply (1) the maximum size of each subsystem is limited, and 
(2) a reservoir is required for both supplemental cooling after water 
is cooled below the specified temperature level and as a place to 
dump reject heat directly from the power plant when no subsystem 
can use it. Reservoir size required for this cooling is shown 
by month in table 2. 
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The LP output gives cost of forcing in activities. This serves 
to rank the relative profitability of enterprises and gives the 
order they would come into the system if more waste heat were 
available. More waste heat woid.d become available if the sizes 
of subsystems currently included in the system were reduced. 

Table 4 [^from ll] shows the greenhouse enterprise is less desirable 
than other possible enterprises including more units of soilwajcming 
and fish culture operated at **non-optimal " temperatures. Greenhouses 
are undesirable because the capital and annual costs given in 
table 5 axe higher than that for other systems, while their rate 
of heat dissipation occupies a middle position. 


INSTITUTIONAL ARRANGEMENTS 

Institutional aspects involve the contractural arrangements for the 
personnel to manage the waste heat system and for the capital to 
constmact and operate the system. Possible options are listed in 
table 10. Costs incurred by a body seeking fee simple acquisition 
include payment of interest suid principal on bonds raised to finance 
ptirchase of the land, administrative costs, cost of compensating the 
affected commxanities for property taxes foregone where land is 
purchased by a tax exempt body aoid leased back for agriculture- 
aquaciilture . Should the utility company decide not to fully control 
the total integrated system and not raise the capital for one or all 
the separate subsystems, it can enter a contractual agreement with 
private entrepreneurs to supply waste heat water. 

The options are limited by two positions held by the electric utility 
industry (1) reliability in the operation of the heat dissipation 

system is essential because the cost of plant shutdown or operation 
at reduced capacity is significant, and (2) direct utility manaigement 
of the agricultural system is neither a goal nor outcome regardless 
of the financial incentives. This viewpoint limits the options to 
either piirchase and leaseback or contractual firm. In purchase and 
leaseback the utility provides the capital and leases the facility 
to a single management firm which maintains and operates all 
subsystems over the lifetime of the power plant. In the contractual 
agreement a cooperative group of firms agree to construct and operate 
the agricultural system with the utility assiming the cost for the 
water transport system. 

Table 6 [[from ll] shows the net present value (NPV) comparisons for 
these three options. In purchase and leaseback free, agribusiness 
gets all the revenue and peiys mansigerial and operating capital costs 
while the utility bears the capital cost for the subsystems and 
water transport system. In purchase and leaseback rent, agribusiness 
pays the managerial costs and costs for operating capital while the 
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utility uses the remaining revenue to defray costs for capital 
for the subsystems and water transport system. In the contractual 
firms option, agribusiness bears all initial capital costs for the 
subsystems and the utility agrees to supply free waste heat eind 
pays the firsm for their dissipation of the waste heat. A payment 
is necessary becavise the difference between the utility operation of 
only the water transport system and the utility operation of the 
total system is negative; hence, a payment would be set by bargaining 
between the utility eind the firm. The details on the methods used 
to calculate the values in table 6 are given in Meekhof ; 

they involve an economic breaJtdown of the terms in equation (1 ) . 


POLICY CONSIDERATIONS 

The MSU study wan limited to heat dissipation, economic, and 
instltutionai arrangements. However, other factors affect the 
desirability of adoption of agricultural usage of waste heat. 

Some of these are contained in a study on energy parks in Michigan. 
Others are shown in studies by others on once throu^ cooling to 
the Great Lakes. 

Economic Expansion 

Beneficial agricultural uses of waste heat were included as one 
collocated industry in an energy park study funded the Federal 
Energy Administration Unlike concurrent energy park studies 

where the park would be located on federal land, the -Michigan study 
assvuned the hypothetical park could obtain privately held land 
near Harbor Beach in Michigan's thumb for the 95^0 acre 23125 MWe 
complex (13 nuclear and 8 fossil fuel units). 

This peripheral involvement in the energy park study occurred relatively 
early in the MSU study. The components suggested and their sizes 
are shown in tables 7 and 8. The sizes were determined by marketing 
constraints. The Initial capital and annual costs were estimated from 
costs for enterprises not \ising wante heat. The heat dissipated was 
estimated from the optimal temperatures for the subsystems. Three 
interesting points emerged from the study. First, over seventy-five 
percent of the waste heat could not be used beneficially either at 
high tempera tiires for industrial and urban uses or at low temperatures 
for agricultiaral uses. Second, the measures for the benefits are 
related to economic expansion either as jobs created or as production 
from energy dependent industry. And third, it was postulated that 
some environmental groups woiold say "... (the park) shoiald be 
accomplished only if the economic and sociad benefits overwhelm the 
economic, social, and environmental costs. The airguments miist be 
^questionably convincing and the documentation must be of unquestion- 
able integrity" . ^ 
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This last point has implications even for the aigricultural system, alone. 
The primary impact in the local economic environment of the atgriciiltural 
system would be caused by the manner in which the land is acquired. 

The conversion of of 9^0 acres of land to purposes other tham it has 
been traditionally used will affect the distribution of wealth and 
Income in the community, the revenue base of the local governmental 
unit, existing input suppliers growth or decay, the long range planning 
incentives of the remaining agricultural land holders, and the existing 
marketing channels for agricultural and aquacultural commodities. 

Once Through Cooling 

The comparisons in table 1 assume the alternative conventional cooling 
is limited to cooling towers. If once through cooling were possible 
several of the results would be changed. The least negative conven- 
tional cooling system woiald be once thro\:igh cooling. Because the 
reservoir would be eliminated, the NPV for the integrated agricultioral 
system would Improve. However, the disch^ge water temperature might 
vary by about fifteen degrees Fahrenheit L^4,Fig.5] which would 
made it more difficult to design the agricultural subsystem to maintain 
aji optimal temperature. 

Environmental groups prevented the operation of the ?00 MWe Palisades 
Nuclear Plant in 1971 imtil a commitment was made to install a 
mechanical draft cooling tower by 197 ^ The setting of agreed 

on environmental standards is a complex issue. Current concern 
includes Cl5il^l*« 

1. a lack of quantitative estimate of potential effects, 

2. dn increase of attached filamentous algae along the 
shoreline, 

3. and the need to assure protection and propagation of 
a balanced indigeneous population of fish and ther 
aquatic species. 

Legislationand court action over the last several years on these 
issues from a national perspective is outlined in a recent Congressional 
Reseeirch Service review C^7J- A lower standeod for thermal discharge 
but not discharge from the agricultural system is permitted if the 
native ecosystem survives intact Reseeirch to quantify the effects 

on once throxigh cooling are in progress. An Argonne study found no 
adverse effects on the migratory behavior of fish tagged at the Point 
Beach Plant [iS]. A continuing MSU Institute of Water Research study 
has found effects imparted to fish, if any, at the Monroe Plant on 
Lake Erie were soon dispersed into the lake population and becsime 
unidentifiable 9] • 
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C0NC3LUSI0NS FOR THE INTEGRATED SYSTEM 

The mathematical varialiLes used contained useful information, jiUttiovigh 
their munerical values are only approximate, Ihiis, the net present 
value criterion is a useful way to structure a study of waste heat. 

The comparative NPV*s do not either clearly prove or disprove the 
concept of waste heat utilization. Because the net present value 
for the integrated aLgricultxiral system is not more negative than 
conventional cooling methods, beneficial \ase of waste heat is not 
clearly disproved. 

Feasible institutional arrangements can be made, but the negative 
NPV implies the electric utility must indirectly pay the agricultural 
comidex to dissipate heat. 

Because water transport NPV dominates the analysis, the primary 
technical problem is the system for water transport, Ihe water 
from the power plant to the 160 acres of fish ponds is carried 
by a concrete pipe 600 feet long and 10 feet in diameter. 

NPV is not a totally acceptable criterion. The NPV certainly 
represents waste heat dependent economic expansion in a new 
industry (catfish culture in a cold climate), or in retaining 
an existing industry (heated greenhouses in the presence of rising 
natural gas prices aind supplies) , but in soilwarming can be 
positive even when no waste heat is utilized, 

A small demonstration project would serve to replace the secondary 
and synthetic data used with more accurate primary data. It would 
also Improve liaison with agribusiness. 


SIMULATION DEVELOIMENT 

The MSU stiidy manipulated mathematicail models rather than executing 
field experiments. An effort wais made to structure the simulation 
and optimization in a manner familiar to the Electric utility 
industry. Desalination of sea water by nuclear power plants ^20] 
seemed to be a problem analogous to the MSU study on waste heat 
utilization. Ibis view was expressed by the MSU group several 
years ago [21]. The einalogy proved to be poor. It is worthwhile 
to outline this formulation of the waste heat problem. Although the 
computer code developed for it was not effective for optimization, 
the code is qirLte useful as a matrix generator for a lineeor (separable) 
programming formulation of the problem. 
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Initial Formiilation 


As an optimization problem the goal xs to maximize the net present 
value. To accomplish this the four type of variables in the NPV 
equation (1) axe manipixLated : R, K, C, and A. The revenue (R) is 
the price of the product time the quantity produced. The production 
in each subsystem depends on its temperature : thus , models are needed 
to predict temperature for a specified quantity of reject heat and 
to predict biological growth rate at a given temperature. At the 
time it seemed reasonable to assume that the agricultural waiste 
heat complex would contain as many as eight subsystems; soilwaxming, 
fish culture, greenhouses, waste treatment, algae culture. Irrigation, 
recreational lake, and animal shelters. Thus, the problem was one 
of allocating the heat between the competing subsystems. It also 
seemed reasonable to assume the initial capital (k) and combined 
annual pumping and operating costs (C) could be found after the 
system was optimized. Also the water transport NPV (Z*) was not 
assignable to any particular subsystem evenubased on the quantity 
of water used. 

To develop the ‘methodology or "tool", soil warming and fish culture 
were chosen as good candidate subsystems. To satisfy the 1983 
zero discharge requirements a cooling reservoir (lake) was included. 
Thus given a power plant with a known quantity of reject heat at 
a known temperature, the problem was (1) to allocate the heat to 
two competing agricultural subsystems on a monthly basis, and (2) 
to determine the size of each agricultural subsystem. Temperatures 
for optimal growth were avaiable in the literature. Applying this 
methodology to all eight candidate subsystems, the optimization 
would give (1) the best size for each subsystem (including some 
of size zero), and (2) would determine whether a particular sub- 
system should be operated at a "non-optimal" temperature in some 
months in order to maJce the best total system. Thus, the best 
uses would be selected and the problem of designing for utilization 
versus dissipation posed by Shapiro would be solved. 

Next ar optimization technique was selected. The mathematical 
optimization technique would need to handle constraints because 
the sume of the reject heat allocated to the subsystems wets limited 
by the reject heat available from the power plant. Because the 
equation to be developed for the heat transfer and biological 
growth had not been decided on a general mathematical programming 
procedure was selected. The "complex" method [.22] met both these 
requirements. The method handles constrained, multi- variable 
problems and only requires computation of functional values, not 
gradients . 
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Subsystems 

The major research effort in small scale waste heat field studies 
is defining the effect of temperature on growth, while a major 
economic concern is to have marketable products from the heated 
enterprises. The same concerns were present diuring the development 
of mathematical models for the revenue (R) in eq^iaation (1 ) . The 
variable R is the first of four types of variables to optimize in 
equation (1 ) . Selecting enterprises for the subsystems was not 
difficult; for candidate "crops" lists of optimal temperature for 
growth wei^ developed and lists of comparative economic advantages 
of production were developed. The crops were the best from the 
combined list. Catfish was chosen for the fish enterprise because 
a detailed economic analysis was available [24]. Tomatoes and 
several other field crop varieties were selected for the soil warming 
enterprise and the economics developed at MSU. Ornamental flowers 
were selected for the greenhouse enterprise and the economics also 
modelled at MSU. 

The price assumptions were: per pound channel catfish $0.30, per 
bushel tomatoes $5.50, ornamental flower rotation and price; 

6-inch Chiysanthemums $2.60, 6-inch Poinsettias $3.25, 6- inch Lilies 
$2.60, 4-inch Geranitims $0.65. 


The heat dissipation was modelled by algebraic energy balance equations 
The optimization method assigned a quantity of heat to be dissipated 
to these equations and the corresponding subsystem temperature computed 
For the fish pond and cooling reservoir, the three equations (I0,15il7) 
in Edinger et al. [25] were solved simultaneously to -yield the thermal 
exchange coefficient and equilibrium temperature. Values for Meyer's 
evaporation term were C.=ll and f =0.00682 + 0.000682 W, where W is 
windspeed in mph. For the soil warming system, the climatic data 
based Kendrick- and-Havens heating pipe model developed by Dewalle [26] 
was used. The heat required to maintain greenhouses at specified 
levels of temperature was computed with the model by Walker [2?] 
with the ventilation lutes adjusted for those used in waste heat 
greenhouses [28] . 

Having established the temperatures, the quantity of crops from the 
subsystems could be predicted. The growth of catfish lased a model 
of the form: 


d W 
d t 


U W 


( 3 ) 


where W is the weight of harvested fish and p. is a parameter modelled 
by a Lagrangian interpolation polynomial using data from Andrews [29] . 
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The growth of field crops used a model from Paltridge and Denhold [,30l 
of the form; 


d G 
d t 

= a G H(s-t) - p G 

( 4 ) 

d W 
d t 

= P G H(t-s) 

( 5 ) 


where G is vegetative material, W is harvested material (tomatoes), 

H is a step function, s is a switch time, and P is related to 
soil temperatiare . Both equations (3) and (4) are reasonable 
approximations to experimental data. Conventional production values 
were used for the greenhouse enterprise. The revenue (R) from 
each subsystem is simply the quantity of harvested material times 
its price. 

Initial Optimization 

The difficulties with the assumptions in the initial formulation became 
obvious when optimization was tried. These were related to: (1) 
the definition of a harvestable crop, (2) design of the water transport 
and (3) institutional limitations on the number of subsystems. 

The formulated optimization problem allocates heat to the subsystems 
on a monthly basis. Only monthly climatic data was readily available. 
However, the harvestable material is only defined in a particular 
month for tomatoes and in a particular size interval for fish. This 
problem was handled by (1 ) optimizing on a yearly rather than a 
monthly basis and (2) using a sensitivity study to define input 
conditions for fish size at harvest. The monthly heat allocation 
pattern was represented as the sum of certain Walsh fimcitions. Each 
function is a pattern with either unity or zero in a given month; 
i. e., to allocate the same unknown heat dissipation to the fish pond 
in the three months Janiiary, February and December the function has 
the form (1,1, 0,0, 0,0, 0,0, 0,0, 0,1). The optimization routine found 
the best amplitudes for these functions. A total of five functions 
was adequate for the allocation to two systems over twelve months; 
thus, the number of optimization variables was reduced from 24 to 5* 

The use of Walsh functions would provide a reasonable approximation 
to a specific optimal temperature pattern obtained analytically 
by others for lobster culture 

Associated with each monthly heat dissipation is a water flow rate 
to each subsystem. When these monthly flow rates were vised to 
determine a water transport system, a more serious problem arose. 

The flow rates are shown in table 9 for three Walsh functions for 
fish ponds auid two for soilwarming. 
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The water transport system detenrmines the next three optimized 
varlahLes in the NPV equation: C, K, and Z*. The ann\oal costs (C) 
are determined hy the flow rate and the friction head loss in an 
established, well-known manner. The flow rates are given and the 
headloss is computed by assuming the pipe diameter and (1) assigning 
a roughness length, (2) using a Moody diagram to compute the 
dimensionless friction factor, and (3) using the Darcy-Weishach 
equation to compute the head loss. This procedure was repeated 
trying different diameters until a satisfactory engineering design 
was obtained. The wide variation in flow rates shown in table 9 
made this design not totally satisfactory. A computer program 
for the calculations involved was found too late in the project 
to be used C32] . Although more designs could have been tried 
using the program, the problem of the variations would have made 
it diffic;alt . to choose correct pipe diameters and piunp sizes. 

In the initial formulation it was assumed that the water transport 
design woifLd not be a problem. In reality the flow rates shown 
in table 9 vary by a factor of 2 for soilwarming and by a factor 
of 5 to 15 for fish pond depending on whether the system operates 
in winter or not. Thus, the water transport costs need to be related 
to the optimized heat dissipation patterns in order to properly 
optimize equation (1). The best way to optimize might be to use 
LP to select from a set of discrete designs. 

The value of the NPV given in equation (1 ) depends on the sizes 
specified for the subsystems (a) which is the last optimized 
variable. Heat dissipation considerations only place an upper bound 
on each subsystem size. The determination of the’ mix of sizes that 
optimize the NPV is determined largely by institutional factors. 

The terms in the NPV equation are basically a subsystem size times 
revenue generated (or indirectly heat dissipated). Some subsystems 
generate large revenues but dissipate heat poorly, others generate 
small revenues but dissipate heat effectively. Thus while revenue 
and heat dissipation considerations could indicate the best mix is 
all greenhouses or all fish culture doing so would possibly mean 
excessively large amoimts of land and/or capital, or large seasonal 
excess capacity in operating capital, labor, or manaigerial skill. 

Also it could create an inability of the marketing structuis in the 
region to handle all the agricultural and aquacultural commodies. 

Th\is, optimizing the NPV eqiaation is not a problem of finding the best 
heat dissipation pattern and compatible size mix, but one of finding 
sizes that are feasible for the institutional arrangement. Table 10 
is a list of possible institutional arrangements. As indicated 
previously the best options are ; purchase and leaseback and contractxial 
rent. 
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SIMULATION STUDY 

The results presented in this paper were obtained using the following 
procedure. The problem of selecting optimal uses was considered 
more as a sensitivity problem than as an optimization problem. An 
LP prograjn was used for the major analysis with technical coefficients 
set from the "complex” program used as a matrix generator. The 
constraint section of the "complex" prograjn was used to set pairs 
of subsystem size and temperature levels. The constraints on heat 
and flow rates were checked with the subsystem models. Associated 
with the temperature levels are crop growth rates and hence revenue. 
Enterprise costs were calculated in detail at a few sizes and 
separable programming was used to interpolate between these values. 


Separable Programmins 


Several related procedures are called separable prograjnming. The 
MSU study used the procedure introduced "by Chames aind Lemke L33D* 
The separable technique was used to (1 ) handle the nonlinearity 
in the cost coefficients and (2) to determine whether subsystems 
should be operated at non-optimal temperature. Other aspects of 
the LP technique were ( 3 ) selecting a reservoir size from the twelve 
monthly values, and (4) combining capiteil amd net revenue internally 
to obtain the net present value. 


The form of the tableau is illustrated in table 11 for a simple 
set of conditions. The system represents three components: a fish 
enterprise operated at the optimal temperature, a fish enterprise 
operates at a non-optimal temperature, and a cooling lake. The year 
has two time periods: summer and winter. The fish enterprises only 
operate in summer. The variables in the tableau are: 

d = Discount factor obtained by summing the opportunity 
cost over the planning horizon, with (R-C) outside 
the summation, 

H = Heat dissipated for a size S in the agricultural 

subsyste, million BTU per hour, 

Q = The quantity of reject heat, million BTU per hour, 

N = The positive net revenue (R-C) from equation (1 ) 

K = Initial capital from equation (1), 

S = Size of the subsystem that detailed costs were compiled, 
it can be zero, 

S* = Size at a level 2 design temperature, 

h = Heat dissipated by the reservoir per acre, million BTU 

per hour per acre. 

The way the four computational features are achieved can be understood 
by studying the tableau. 
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Severai aspects of the tableau and programming method are of interest. 
The heat dissipation rates, H, are non-negative; that is, heat is 
transferred from the subsystems to the environment, not vice versa. 

Hence the maximm size of each subsystem is limited by heat dissipation 
considerations since the quantity of heat rejected from the power 
plant is limited. More importantly the agricultural subsystems only 
contribute to the net present value if they also dissipate heat. 

The size of the agricultural subsystem is limited to be between 
sizes and by the separable constraint; thus, should be set 
to zero if the '^subsystem if not necessarily in the integrated system. 
Unlike the "complex" program where the subsystem costs coiild not be 
related to the water transport costs and design, different discrete 
water transport designs can be specified as level 1 and level 2 
subsystems in the IP program. Ihe "complex" program has some feaiures 
that can be used to set technical coefficients: given a quantity of 
heat to be dissipated it will check the implicit constraints on 
power plant temperature and on designed pimping rates and obtain 
a feasible subsystem operating temperature for heat transfer consider- 
ations. 

Program Notes 

The mathematical methods developed by the HSU group and those developed 
by the Price group at Comell[3^1 have helped to isolate the key 
variables in the waste heat problem. 

Several aspect and potentials of the LP programming method beyond 
the goals of the MSU study should be noted. 

The technique used for internal discounting has been used in another 
context to weight monetary and pollution costs 1_35D* This could be 
used in waste heat by replacing the variable S with a measure of 
production of catfish wastes and a weight analogous to the discount 
put in the objective function in order to find a least cost, least 
pollution integrated system. 

The pumping costs imply that the benefits from the agricultural system 
is derived at the cost of using limited energy resources to pump the 
water. This treide-off has been discussed in a nice mathematical way 
in the context of altering the | thermal plume by pumping by others • 
In the LP analysis the variable S could be replaced by a measure of 
energy usage and the weighted NPV obtained. 

Finally, linear (separable) programming is also a method familiar to 
the electric utility industry in the context of designing a dry 
cooling tower and better for the waste heat problem than [20,21]. 
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TABLE 1 

COMPARISON OP COOLING METHODS 
NET PRESENT VALUE (Z) 

(IN THOUSANDS OF DOLLARS) 


Discount Rate 
and 

Planning Horizon 

Agri cul tttral 
Method 
at Optimum 

Wet 

Mechanical 

Draft 

Wet 

NatureJ. 

Draft 

10 ?S at 
30 years 

-43395 

-56912 

-64344 

12 ^ at 
30 years 

-38783 

-52123 

-59255 

15^ at 
30 years 

-33708 

-46923 

-53730 

lO^ at 
25 years 

-42255 

-55690 

-63046 

12 ^ at 
25 years 

-38092 

-51382 

-58423 

15% at 
25 years 

-33465 

-46568 

-53352 


Month 

TABLE 2 

MONTHLY RESERVOIR REQUIRED 
FOR SUPPIEMENTAL COOLING 
(ACRES) 

Area Month 

Area 

January 

151 

July 

350 

February 

143 

August 

357 

March 

195 

September 

357 

April 

200 

October 

242 

May 

313 

November 

218 

June 

329 

December 

154 
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TABLE 3 

SHADOW PRICES FOR OPTIMAL SYSTEM 
(IN THOUSAND OF DOLLARS PER ACRE) 


Fish Pond Activity 


10% at 
25 years 

96.6 


15 ^ at 
25 years 

69.2 


10 % at 
30 yeaJTS 

100 . 


15 % at 
30 years 

69.9 


Soilwarming Activity 1616. 


937. 1713. 


952 . 


Fish Pond Resource 
Soilwarming Resource 


39.3 

17.1 


27.3 

11.1 


40.9 

18.6 


27.7 

11.3 


(IN THOUSAND OF DOLLARS PER MILLION BTU) 


September Waste Heat -0.975 -O .766 -1.01 -0.774 

Resource 


TABLE 4 

COMPETITIVE POSITION OF NON-BASIS ACTIVITIES 
(IN THOUSANDS 0F DOLLARS PER ACRE) 


Order 

Activity 

Cost 

5 

Tomatoes, 75 Acres at 
10 % Below Optimal 

282 

6 

Tomatoes, 75 Acres at 
Optimal Temperature 

295 

9 

Fish Ponds, I 60 Acres at 
10% Above Optimal 

482 

14 

Greenhouse, 1.25 Acres at 
Optimal Temperature 

622 
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Subsystem 


Fish Pond 
Soilwaxming 
Greenhouse 
Reservoir 


Rate 

and 

Horizon 

10 ^ at 
30 years 

12 ^ at 
30 years 

15 % at 
30 years 

10 % at 
25 years 

12 % at 
25 years 

15% at 
25 years 


TABLE 5 

ECONOMIC COEFFICIENTS IN NPV EQUATION 
FOR OPTIMAL SYSTEM 
(IN THOUSANDS OF DOLLARS PER ACRE) 


Initial 
Capital , K 

11.9 

7.75 

260. 

10.5 


Anm:al 
Costs, C 

11.4 

2.09 

239. 

1.61 


Annual 
Revenue , R 

16.5 

4.18 


ZERO 


TABLE 6 

MANAGEMENT COMPARISONS' 
NET PRESENT VALUE 
(IN THOUSANDS OF DOLLARS) 


Total 

Utility 

Operation 

Purchase 

Leaseback 

Free 

Purchase 

Leaseback 

Rent 

Utility 

Water 

Transport 

- 4339*5 

-60542 

-53146 

-39967 

-38783 

-53372 

-47052 

-34905 

-337O8 

-45623 

-40471 

-29424 

-42255 

-60738 

-53616 

-40676 

-38092 

-52293 

-46139 

-34123 

-33465 

-45090 

-40018 

-29034 
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TABLE 7 

COLLOCATED INDUSTRY AND AGRICULTURE 
HYPOTHETICAL ENERGY PARK AT HARBOR ^CH 


Industry 


Capacity 


Paper Mill 
Chemical Complex 
Petrolevim 

Coal Gasification Plant 
Mini Steel Mill 


2.000 tons per day 

$240 million annual sales 

250.000 hhls per day 

250 million cfd of gas (300 Btu/cu ft) 

150.000 tons per year 


Biocomplex 
Irrigation 
Catfish Cultiire 
Greenhouse 

Grain Drying, Low Temperature 
Waste Treatment 


Capacity 

3200 acres outside park 
500 acres in 20 acre ponds 
300 acres (vegetahles and flowers) 
20 million bushels per year 
85 acres of algae ponds 


Park City 

District Heating 135»000 people in lO square miles. 


TABLE 8 

PARK REQUIREMENTS 

HYPOTHETICAL ENERGY PARK AT HARBOR BEACH 


Measure 

POWER 

PLANT 

INDUSTRIAL 

PARK 

BIO- 

COMPLEX 

PARK 

CITY 

TOTAL 

Land 

(acres) 

/95IO 

3000 

960 

6400 

19870 

Capital 

(million dollars 

1975) 

, 10490 

1000 

59 

• • 

11549 

Operating Labor 
(man equivalent) 

1600 

6400 

300 

• • 

8300 

Water Usage 
(million gpd) 

170 

52 

27 

• • 

249 

Heat Dissipated 
(million BTU/hr) 

11 7800 

14000 

19200 

2850 

153850 


Source; [13] volume 1 , tables 1 and 24 . 
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TABLE 9 

MONTHLY FLOW RATES FOR OPTIMAL SYSTEM 
(IN THOUSANDS OF GALLOl^ PER MINUTE) 


Month 

Fish Pond 

Soilwarming 

Cooling 

Reservoir 

January 

463.6 

• • • 

86.4 

February 

342.1 

• • • 

207.9 

March 

130.0 

• • • 

420.0 

April 

128.9 

• • • 

421.1 

May 

29.9 

10.9 

509.2 

June 

31.1 

6.2 

512.7 

July 

32.6 

6.3 

511.1 

August 

33.2 

• • • 

516.8 

September 

33.9 

• • • 

516.1 

October 

271.6 

• • • 

278.4 

November 

170.9 

• • • 

379.1 

December 

527.8 

• • • 

22.2 


TABLE 10 

SITE ACQUISITION AND MANAGEMENT OPTIONS 

Fee Simple Acquisition 

I^archase and Manage 

Purchase and Leaseback 

IHirchase and Resale on Condition 

Less than Fee Simple Acquisition 
litrchase Easements 

Contractural Agreements - For Rent liroperty Interest 
Waste Heat Water Cooperative 
Contractual Arrangement 

I^iblic Authority 
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TABLE 11 

SOTEMATIC BEPRESENTATION OF 


Objective Function 
Winter Reject Heat 
Summer Reject Heat 
Winter LaJce Heat 
Summer Lake Heat 
Capital Transfer 
Revenue Transfer 
Lake Cost Transfer 
Level 1 Separable 
Size Transfer 
Level Z Separable 
Size Transfer 
Separable Mix 


0 0 0 0-1 d-d 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

1 0 0 0 0 0 0 

0 1 0 0 0 0 0 

0 0 0 0 -1 0 0 

0 0 0 0 0 -1 0 

0 0 0 0 0 0 -1 

0 0 0 0 0 0 0 

0 0 -1 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 -1 0 0 0 

0 0 0 0 0 0 0 


THE SIMPLEX TABLEAU 

0 0 0 0 0 0 0 

1 0 0 0 0 0 0 

0 1 0 H^ HJ- H| 

1 0 -hj_ 0 0 0 0 

0 1 -h2 0 0 0 0 

0 0 k^ K, Kf K* 

0 0 0 N* N* 

0 0 C 0 0 0 0 

0 0 0 1 1 0 0 

0 0 0 S^ S^ 0 0 

0 0 0 0 0 1 1 

0 0 0 0 0 S|^ S| 

0 0 0 1 1 1 1 
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